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Non- Fermi liquids can be studied using holographic duality. The low energy physics of 
a holographic Fermi surface is controlled by an emergent scale invariance. After reviewing 
these developments, we generalize the holographic calculation to include in the bulk 
action the leading irrelevant operator, which is a dipole coupling between the spinor 
field and the background gauge field. We find that this dipole coupling changes the 
attainable low-energy scaling dimensions, and changes the locations of the Fermi surfaces 
in momentum space. The structure of the holographic framework for non-Fermi liquids 
is, however, robust under this deformation. 



O^l 1. Introduction 

> 

Holographic duality is a wonderful development arising from string theory which 
offers a new perspective on strongly coupled quantum systems. The duality solves 
certain such systems in terms of an auxiliary theory of gravity in one extra dimen- 
sion. This extra dimension plays the role of the renormalization group scale, and 
Einstein's equation is an RG evolution equation^] 
y—i The question "To which strongly coupled systems does the duality apply?" has 

not been settled. The best understanding at the moment is for certain relativis- 
tic, supersymmetric gauge theories. Though these systems are not yet known to be 
realized in Nature, this shortcoming can be mitigated by judiciously chosen ques- 
tions about long-distance physics, where the same universal behavior can arise from 
systems which are very different microscopically. 

The basic dictionary between the 'boundary' quantum system and the 'bulk' 
gravity theory is as follows. The gravity theory is classical when the boundary 
system has a large number (which we will call N 2 ) of degrees of freedom per spatial 
point. There is a one-to-one correspondence between local operators in the boundary 
theory and fields in the bulk. Correlation functions of these local operators may be 
computed by solving classical wave equations for the corresponding bulk field. To 
study thermodynamic equilibrium of the boundary system, one places a static black 
hole in the bulk geometry. 



X 



a For reviews of these developments in the spirit of the present work, see fljo 
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Attempts have been made to apply this technology, with some success, to the 
quark-gluon plasma, quantum-critical transportpl] and ultracold atoms at unitar- 

ityjaisa 

One is led to wonder whether this apparatus can be applied to the study of 
metallic states. Most of the present understanding of metals relies on Landau's 
Fermi Liquid ansatz: one supposes that the low-excitations are long-lived electron- 
like quasiparticles near the Fermi surface in momentum space. This ansatz results 
in a powerful and robust low-energy effective field theory with many successes, but 
there are systems to which it does not apply. The feature of such "non- Fermi liquids" 
on which we focus is the presence of a sharp Fermi surface of gapless excitations, 
which lack a description in terms of long-lived quasiparticles. 



In 11 - 15 , it was shown that holographic duality can describe Fermi surfaces. 
Indeed, they arise from the minimal ingredients necessary to pose the problem. 
The system under study is, microscopically, a 2+1-dimensional relativistic confor- 
mal field theory (CFT); this means that the dual geometry is asymptotically four- 
dimensional anti-de Sitter space, AdS^. Assuming that this CFT has a conserved 
U (1) current (a proxy for fermion number), the gravity theory must include a mass- 
less photon field A„. In order to study fermion response, the CFT must contain a 
fermionic operator, with some charge q under the U(l) current, and some scaling 
dimension A at the short-distance (UV) fixed point. To introduce a finite density 
of U(l) charge, one studies the charged black hole in AdS^. This finite density 
breaks the Lorentz symmetry and scaling symmetry of the boundary theory. The 
zero temperature groundstate is described by an extremal black hole. 



In Ref. 14 it was shown that this finite-density ground state exhibits an emer- 
gent scaling symmetry which is manifested by the fact that the near-horizon region 
of the extremal black hole geometry is AdS-z x M. 2 . By holographic duality, this 
geometry is dual to a fixed point field theory, which we call the IR CFT. The IR 
scale transformations act on time but not space, and in this sense the IR CFT has 
dynamical exponent z — oo. This emergent quantum critical behavior governs the 
low-frequency fermion response. The retarded fermion Green's function exhibits a 
Fermi surface, near which it takes the form 



W(kM-k)) = , hl t7 2v (l) 

k i — — u) — hoe^uj^ 



with kj_ = \k\ — kf, and real constants hi, vt , 7, v. v is a scaling dimension in the IR 
CFT. In this way, the emergent conformal symmetry controls the dispersion of these 
excitations. This duality construction provides a large theoretical playground of non- 
Fermi liquid fixed points, giving a handle on a difficult strong coupling problem. In 
particular, the case with v = 1/2 has several features in common with the "strange 
metal" phase of high T c superconductors. 
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In this workup | we consider effects on these holographic Fermi surfaces of an 
additional dimension five operator 

$(g m + g e T)T MN ^F MN (2) 

in the bulk action corresponding to magnetic and electric dipole moments for the 
bulk fermions (T = T^T^Y*, T MN = \ [r M , T N }). Our motivation for this addition 
is several- fold. First, we would like to investigate the robustness of the previous dis- 
coveries. The choice of action for the bulk fields in the calculations described above 
was motivated by Landau-Ginzburg- Wilson Naturalness: the lowest-dimension op- 
erators respecting the symmetries were used. This Naturalness criterion is usually 
enforced by the renormalization group. In the bulk gravity theory, the status of 
this principle is not clear because of our poor understanding of quantum gravity. 
This motivates an exploration of the sensitivity of the results to RG-irrelevant bulk 
couplings. 

In particular, our previous calculations used the canonical Dirac action. The two- 
point function is insensitive to self-interactions of the spinor field at leading order in 
the 1/N expansion. However, there exists a large class of higher-dimension operators 
which are quadratic in the spinor that can possibly change these conclusions. This 
work on dipole couplings is an attempt to investigate systematically the effects of 
such operators. Do these higher dimension operators in the bulk drastically alter 
the existence of Fermi surfaces in the boundary? We will find that they do not, but 
rather the main effects of the dipole couplings are to change the IR AdS2 scaling 
dimensions and to change the locations of Fermi surfaces in k space (which we will 
find numerically^ As such, we have constructed a much larger parameter space of 
non-Fermi liquids for study. 

The dipole moment couplings are a natural starting point for an exploration 
of irrelevant operators, since the structure of the calculation is largely preserved. 
Further, they are generic in the following sense. The AdSi under study frequently 
arises as a factor in a higher-dimensional bulk spacetime AdS$ x X, where X is 
compact. There will be an effective theory of the light modes (or a "consistent 
truncation") on AdS±. Even when one starts with the simplest Dirac Lagrangian on 
AdS\ x X, if the charge on the black hole comes from the Kaluza-Klein gauge field, 
such dipole terms generically appear in the low energy effective theory on AdS± Q 
These couplings also appear in some string theory realizations of the ingredients 



b Note Added 1: After this paper was completed, but several months before it appeared on the 
arXiv, the papers [44|45] appeared, which are the first published papers to study the effects of the 
magnetic dipole coupling g m . We feel that our perspective on the subject is still worth sharing. 

I do not study Re- 
tention to 2+1-dimensional field theories, in which case the gravity theory lives 



The papers |44|45| do not study g e 
c We restrict our att 



in four dimensions. 



d 



Note Added 2: However, 



44]45 arger values of p m than are 

studied here, Fermi surfaces are not present, and further that the low frequency spectral weight is 
suppressed. We make some comments on these points at the end of this paper. 
e We thank N. Iqbal for this point. 
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described aboveP^^ 

2. Fermion Green's Functions from Holography 

In the next three sections, we review previous worlii21iil22l on the holographic compu- 
tation of two-point functions of fermionic operators in a 2 + 1 dimensional boundary 
CFT with a finite U(l) charge density. For simplicity, we work at zero temperature. 

In the bulk, this ensemble corresponds to a black hole in AdS^ charged under a 
U(l) gauge field. We employ the bulk action 



Sbulk = / \'-</ 



R 2 9 F 



(3) 



where R is the AdS radius, k the Newton's constant and g F the gauge coupling. 
More specifically, the relevant solution is the charged AdS± black hole, 

d s 2 = ^(-fdt 2 + dx 2 ) + ^, f = i + &-M A t = - — ) (4) 
R l r A j r 4 r 6 r 

with Q,M the black hole charge and mass respectively and /x = g F Q/ (R 2 Tq). tq 
is the outer horizon, i.e. the largest solution to f(ro) = 0. In the boundary, this 
geometry corresponds to a theory with finite charge density and temperature 

2Q = _3r^ ( _ Q 2 

n 2 R 2 g F ' AttR 2 V irl 

At extremality, the inner and outer horizons merge into a double zero of / and 

M = 4 C^j / , Q = VSrl => T = . (6) 

We will work in units with H = 1. In addition, in our numerical work, we will often 
put the horizon at ro = 1 and set g F = 1. 

We want to study the Dirac equation in the bulk, 

T m DmV - mtf = . (7) 

Here T M is related to the usual flat space gamma matrix by a factor of the vielbein, 
pM _ e A/^pM anc j £) M — q m _|_ 4ci; a hMr ab — iqAM with w a ijif the spin connection. 
One can nicely cancel off the spin connection contributions to this equation by 
defining 

* = {-gg rr )-*e- iut+ikixi ilj (8) 

with x — (x, y) the spatial directions on the boundary. Substituting and rearranging, 
we get 

' 9i 



-(ra r -mV^)^ + iif,r^ = 0, K„ = {-u,ki) (9) 

Qrr 

with 



J^-^ + fX^l-y)). Ii, = l"i- (J.0) 
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This system of four coupled equations becomes simpler by rotating the k momentum 
to be entirely in the x direction (which we can do by rotational invariance) and by 
a choice of gamma matrices adapted to this frame, 



-a 3 

-cr 3 



ia L 
ia 1 



-a 2 
a 2 



-ia 2 
2 



ia 



(11) 



Defining ip 



$ 2 



and rearranging gives 




ia z <f> n + J— k{-l) a a l <$> a (12) 
9v. 



with a = 1, 2. This gives two decoupled, real 2x2 equations 
like 



A solution of the Dirac equation, in the basis ( 11 ), behaves near the boundary 



$ Q - a a r m ( 1 ] + b a r" 



(13) 



In terms of the eigenspinors of r r , tp + , tjj_ , we have 
V>+ - A(k)r m + ... , A(k)=( ai ); V- ~ D(k)r~ 



These coefficients are related by a matrix S, 



D(k) 



S 



S3 S4 



(14) 



(15) 



Since the two a equations are decoupled, we can choose independent boundary 
conditions that do not mix $i and $2 giving S2 = S3 = 0. The standard prescription 
for calculating the retarded Green's function gives 



G R = iS 7 * = 



with 



7*' 
7*' 



61/01 
62/02 



(16) 



(17) 



where U is the basis change between the $ a basis and the chiral basis. 

For < m < 1/2, there exists an inequivalent holographic prescription using this 
same bulk action which describes a different boundary theoryPS In this "alternative 
quantization," the roles of the source (a) and response (6) are switched, and similar 
reasoning leads to G R = — 1/G^. The boundary CFT in alternative quantization 
flows to the usual one upon adding O^O to the CFT lagrangian, where O is the 
operator dual to if). 
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3. Low Frequency Limit of Retarded Green's Functions 

We are interested in frequencies small compared to the chemical potential, \i. 
Naively, we should expand in a series in lo. However, the term proportional 



to uj in (12 1 is dominant at the horizon and so we cannot treat w as a small per- 
turbation there. To deal with this, split the r-axis into two regions, an inner region 
(with coordinate £) and an outer region (with coordinate r). The inner region is 
specified by: 

r-r Q = , e < C < oo (18) 

and the outer is 

<r-r (19) 

with i?2 = l/v6 (recall that we have set the AdS radius, R = 1). The strategy now is 
to develop the solution as a perturbation series in u> with Q in the inner region and r 
in the outer region. Because the distinction between inner/outer involves ui, the inner 
region equation no longer blows up in the u) — > limit and the perturbation series 
between the two regions is reshuffled. The statement that results are independent 
of the matching point e is a holographic version of the Callan-Symanzik equation. 

Let us examine the lowest order solution in the inner region by taking the limit 
lu — >• 0, e — >• 0, luR\/c — > 0. In this limit, the Dirac equation is 

(-flfc + ^V)$ a = (a + c ^W* a + ^ ( _i r iv$ a (20) 

C C C r 



with e3 = 3f/v 12. This is precisely the Dirac equation for a spinor in AdS?. x R : 

ds^^i-dr' + dC^ + ^d^ (21) 

in the presence of a constant electric field e^. Here iji a = (— gy^) -1 / 4 ^, where £ = 
is the boundary, and the AdS2 time coordinate is r = Xt. To be consistent with the 



definition of £ (18) the parameter A should be set equal to cj; the parameter A is 



introduced to avoid the awkwardness of rescaling the time coordinate by a frequency. 



Relatedly, in (20) we have defined Q = w/A; this is the frequency conjugate to the 
AdS2 time coordinate. 

We are interested in matching to the outer region where £ — > 0. Here, the 
equation becomes 

/ -mR 2 e 3 q+ (-l) ah ^\ , . 

This matrix has eigenvalues ±z/ Q with 



+ m?Rl - q 2 ej . (23) 
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v a is the scaling dimension in the IR CFT of the frequency-space operator dual to 
ip a . The corresponding eigenspinors are 

V ±<* = L „ I l\akR 2 ) • (24) 



xe 3 <?-(-l, ro , 
At the boundary of AdS%, the solution is therefore 

^=v_ a C Va +G^)v +a C (25) 



Generalizing ( 16 ), the AdS 2 Green's function (in the presence of constant E field) is 



Ox in (25). To get the retarded function, we must set infalling boundary conditions 



at the horizon. The AdS 2 Dirac equation (20) is solvable} 2 ^" and the associated 
retarded Green's function is 

R[ ' ' T{2v a )Y{\-v a -iqez) ( m + ±)R 2 - iqe 3 + v a K 1 

(26) 

Note that by momentum conservation in R 2 , operators with different k do not mix. 

Now we look at the outer region equations. Here, we can safely set w = to 
get the lowest order solution; a basis of solutions is specified by the IR boundary 
condition 

v¥i~v±J-^-) ±Va (27) 



r - r 



Matching, we conclude 



*2 (0) =^ + aS(^. (28) 



Now, in the outer region, we can perturbatively expand the linearly independent 
solutions 

V± a = V ( ±l+^V ( il + ■■■ (29) 

where we have already solved for 77^. The higher orders can be obtained by solving 
the Dirac equation and requiring that the solution has no piece proportional to the 
lower order ones. Thus, the matching is entirely determined by the lowest order and 
we conclude 

*%=v+ a + Q%(f>)v- a (30) 

To know r)± a we must solve the Dirac equation everywhere - we must have all the 
UV data. G%{w) is however determined entirely by the IR CFT. At the boundary 



of AdSi, we have (from (13)) 



d B 2~«£2r-(;)+6£1r-(j) (31) 



giving the full Green's function perturbatively using ( 16 1 
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4. Fermi Surfaces 



Let us suppose that there exist certain kf where a?„(fc/) = 0. This will only happen 



for real v a . For small k± — k — kf and small u), the Green's function (32) can be 
written 

G£(w, k) pa ^— = (33) 

-L v f 



with the (real) constants h\.2,Vf in (33) determined by the UV data a±\ and 
7fe, is the phase of Qn(kf,oj). This Green's function has a pole in the complex w 
plane at 



w c = w*(fc) - iT(A) = H^J *' e " R/ _ < 2 . (34) 

[vfk ± - v f h 2 e lk f (vfkj_) 2Uk f v kf > \ 

We interpret the uj = 0, k± = singularity as a Fermi surface and the finite w poles 



as particle-like excitations above this Fermi surface. Looking at ( 34 ) , the excitations 
have dispersion relation (k) cx k^_ and widths T(k) oc k s ^ with 

and 5=]^ . (35) 

For < ^, the width and energy are comparable and the excitations are not 
stable quasi-particles; these are non- Fermi liquids. For v^ F > ^, as we scale towards 
the Fermi surface, the ratio of lifetime to energy goes to zero and these are stable 
particles. For y kp — 1/2, both Q%(u>) and a^ a have poles which cancel, leaving a 
log in the Green's function. It is 

G%(u, k) « ^— (36) 

■ + CiU) + CiwlogW 

with ci real and c\ complex. Such a Green's function is characteristic of "marginal 
Fermi liquids," which provide a phcnomcnological model of the strange metal state 
of the cupratesP^l 

Thus, we have obtained Green's functions for a family of excitations about a 
Fermi surface. The low energy properties, such as the form of the lifetime and 
dispersion relation are entirely determined by the scaling dimensions of an emergent 
conformal field theory. 



5. Turning on Dipole Couplings 

We want to look at the effects of changing the intrinsic electric or magnetic dipole 
moment of the bulk spinor on the structure of these holographic Fermi surfaces. To 
do this, we use the bulk Lagrangian density 

C = i(¥ r M D M V - m^) - *{g m + g e T)T MN ^F MN (37) 
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with (in our basis (11)) 







-1(7 





(38) 



The Dirac equation is now 

(T M D M - m + i(g m + g e r)T MN F MN ) * = (39) 
Once again^we can cancel the spin part of the covariant derivative by making the 
definition i s). Using F rt = [ir /r 2 , and the definition ([8]), with ip = ( ^ ) in the 



basis ( 11 ), we get 



Mg (l-~))^$ a 



-9tt \ r / V gi, 

+ 2^r (~g u y^{gm(T 1 ^ a + g e cr 3 $ (3 )^ r 



(40) 



where again a = 1,2 and the index /3 = 3 — a awkwardly indicates the other 
component. Note that when g e ^ 0, the Dirac equation is no longer block diagonal 
in this basis, though it is still real. The dipole terms have no effect on the boundary 
behavior of this equation. 

However, because there is mixing when g e ^ 0, the process for extracting the 



Green's function is slightly more complicated. Equation (151 still holds, but we 



can no longer choose two sets of boundary conditions such that Gr is diagonal. 
Instead, we use two sets of linearly independent boundary conditions, /, and 7T. 



(15) becomes 



b{ b{ 
bl hi 1 



si s 2 

S3 S4 



(41) 



or B = SA in matrix notation. The Green's function is Gr(lo, k) = S = BA 



5.1. Discrete Symmetries 

We can discover several discrete symmetries by examining the effect of conjugating 



the Dirac equation ( 40 1 and the infalling boundary conditions with certain sim- 



ple matrices, U. For example, when our two sets of infalling boundary conditions 
correspond to the two different spins, conjugating with the matrix 



U = 



1 

1 



(42) 



switches the sign of k in the Dirac equation and switches the two sets of boundary 
conditions. We learn that 

G(oj,-k) = UG(u,k)U (43) 

When g e = 0, we can take a diagonal basis, leading to Gi(uj,~k) ~ G2(w,fc). For 
the general mixed case, we note that det G(u>, —k) = det G(tu, k), so that our graphs 
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of Fermi surfaces in the (k, q) plane will be invariant under k — > —k. In a similar 
way, by examining the effect of U on (40) and on the boundary conditions, the 
choice 



U = ( U q ^3 ) = rn =^ G(-u),-k;-q,-g m ,g e ) = -G*{w,k\q,g m ,g e ) . (44) 



This, along with the first discrete symmetry, implies that our Fermi surface plots 
with g m = will be symmetric under q — » — q. Finally, the choice 



U= ( a ® 2 ) =r-r i =^ G(uj,-k;-m, -g m ,-g e ) = - [G(uj, k;m, g m , ge^' ' 



(45) 



In particular, this implies that switching to alternative quantization is equivalent 
to taking (m,g mi g e ) (-to, -g m , -g e ). 



5.2. The Low Frequency Limit 

Again we develop a perturbation series in uj by splitting the r-axis into inner and 
outer regions. The lowest order inner region equation is 



mR 2 
-d ( + — — cr J )$ c 



ge 3 

c 



ia 2 <P r 



+ 2^(g m <T 1 <S> a +g e cr 3 $0) — 

C it2 



fc(-l)<V$ c 



f 



Near the boundary of AdS2 , we get 



(d c <S> a = ia 2 qe 3 $ a - R 2 (toct 3 + A/qO- 1 j $ Q + 2 ffe -g-<7 3 $ /3 



with 

M a = _(-l)»A+2^ ftn . 
Again —Cd^tp = U(g e ,g m )4> and the four eigenvalues of U are 



(46) 



(47) 



(48) 



±- 



1 



\ 



k 2 



k 2 , 



(m 2 + )jR 4 + e 2 (4(5 2 + 3 2J _ q 2 R 2 } ± ^ R 2 g 2 + g 2 (m 2 + ) 



(49) 

where the 1,2 correlates with the ± in the square root. Thus the dimensions of 
operators in the IR CFT are significantly changed. In the case of g e = g m = 0, the 
usual case of two degenerate eigenvalues obtains. 

By making a basis change on ( 46 ) , we can block diagonalize it (though we cannot 

(50) 



do so for the full Dirac equation): 
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with 



\ 



4e§ (g 2 +9 2 m )R 2 + 


I™ 2 + S) 


(g 2 m $ o +9l 








4e\ (g 2 +g? n )R 2 + 


r 2 + S) 


R G 2+^e 2 


[gl^+9 2 e 


(™ ,+ *) 





R-2 



(51) 



This is exactly the same AdS 2 Dirac equation as (20), after the replacement 

k v± 

To R2 

m — > 



(-ir- 



v\ - q 2 e\ = vx,i 



(52) 



Thus, 

g£ a M = e— ^"^ffi 1 + y : ( " : :' x ~ = ^ 2 . ( 53 ) 



r(2n, 2 )r(l - z/i,2 - tgea) 



-zf± - i?e 3 + 1/12 



As in (28 1, we can match in the outer region onto either G R {u>) or ^(w), this defines 
our two boundary conditions. The components in the outer region, however, will 
generically be mixed 4-spinors. We will have two solutions 



We can expand the matrices A and B (41) perturbatively in ui near the boundary. 
For example, 



B = + uB^ + 0{lo 2 ) + (Bi 0) + ujB W + 0{u 2 ))g R {u) 



with 



Qr{uj) ={ g%(w) 

The equation for the low frequency Green's function is (to order u 2 ) 



(55) 
(56) 



Bf +u J B^ ) + (B i ° ) +LoB {1) )g R (u J ) = G R {u, k) [{Af + + (A m + uA^)Q R {u) 

(57) 



All previous equations for correlation functions (such as (33)) hold, with a a 's and 
6 Q 's replaced by matrices A and B, Gr(u)) replaced by the matrix (53), and all 
denominators replaced by matrix inverses. The Fermi surface is now defined by 

det[^ 0) (fe/)] =0 . (58) 



The dispersion relation and width - the analogues of ( 35 ) - are determined by 
solving 



det 



(Af(k f ) + d k Af(k f )k ± +ujA^(k f ) + A W (k f )g R (u) 



. 



(59) 
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Thus, in general, the dispersion relation and width will be controlled by the smallest 
of v\ , v 2 . 

For simplicity, we will find it easiest to deal with nonzero g m and g e separately. 
5.3. g m y£ 0, g e — 

For g e — 0, there is no need to do any basis changing; the Dirac equation is block 
diagonal. The near horizon equation is (setting ro = 1) 

i-Oc + = (a, + ^)ia^ a + + 2 ^)a 1 $ Q . (60) 

Thus the effect of the magnetic dipolc in the near horizon limit is to shift the mo- 
mentum oppositely in the two blocks. In the matching region (the AdS2 boundary), 
the solution goes like 



with 



and 



® a ~ C-V+a + C a V-a (61) 



R> v '(kR 2 2 + (~l)«2 ffm e 3 ) 2 + m^R\ - q^e\R\ (62) 



( ™R 2 ~F v a \ 

v±a = W q -^)-(-irkR 2 ) • (63) 

Matching onto the near horizon region, 

$°=»? + a+aS(w)»?-a (64) 

where the effect of the dipole coupling is to shift k in G%(u), and to change the UV 
data rj± . The AdS2 Green's function is 

G a (oj) = e~ iirVa F (~ 2l/ ") r ( 1 + Va Z iqe ^ x ( TO + ~ ige3 ~ Va (2u) 2va 
r (2i/ Q )r (1 — v a — ige 3 ) (m + in„) — ige 3 + v a 

n a = (-irk+ 2 ^ (65) 
5.4. g e ^ 0, g m = 

In this case, the (ro = 1) near horizon equation can be block-diagonalized to the 
fornJO 

3 



- a C l> Q = (W + ?)^2$a + y((-l)% 2 + ^ + "^) CT l*a • (68) 



f In the basis |^U , the Dirac equation takes the form 

(-9< + ^^(J3)*a = (w + ^r)ia 2 3> a + ^(-l) a ka^ a + ^g^^f) (66) 
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In the matching region, the solution again goes like $ Q ~ £ Ua v+ a + £ Ua V- a with 



v a = ^^W™ 2 + k 2 R\ + (-l)«2 5e e 3 ) 2 - q 2 



(69) 



V± a 



e 3 (9 - ^) - (-l) a Vm 2 + k 2 R 2 



The full outer region solution is $^ = 7] +Q + Q^(u>)r)- a with 

r(-2^)r(i + 1/, 



r(2i/ Q )r(l - f Q - ige 3 



«ge 3 J ^ zn a - iqe 3 - v a 



in a — ige 3 + v a 



{-l) a yjm 2 + k 2 + 



2e 3 g e 
i?2 



(70) 



6. Numerical Results 



To find Fermi surfaces, we look for kf such that a? (kf) — 0. By (31 1, this corre 



sponds to uj — solutions to the Dirac equation which are normalizable (because 
of mixing, the process is slightly more complicated for g e ^ 0; we review it below). 
We implement this procedure by numerically integrating the to = equation to the 
boundary and looking for zeros of cffl for some range of k and q. 

Such numerical work was previously dondHEH for g m = g e = . There, it was 
found that Fermi surfaces existed in branches in the (k, q) plane that were basically 
straight lines jutting out of an oscillatory region (a region where the AdS2 operator 



dimensions ( 23 ) are imaginary and inside which there exist no Fermi surfaces) . See 
Figure [T] for such a graph with m = 0.4. The oscillatory region is shaded green. 



Fig. 1. Fermi Surfaces g e = g m = 0, m = 0.4 



The unitary transformation which block diagonalizes this is <J> = 



U 



2s/2\Jk 2 + m 2 + mVrn 2 + k 2 



(-A+ -A. 
A- -A 4 
-A_ -A-f 
V A+ -A. 



A- -iA+\ 
tA + iA- 
iA + iA_ 
iA- -iA+ J 



A± = m ± k + yrt 



- k 2 



(67) 
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6.1. g m # 0, g e = 



For g„ 



7^ 0, the above structure is preserved; there are Fermi surface branches 
jutting out of oscillatory regions in the {k,q) plane. By (62). turning on g m keeps 



intact the shape of the oscillatory region, moving it to larger k (it also moves 
another copy associated with the opposite spin to smaller k, but we focus on k > 
as everything is k — ¥ —k invariant). We make the following qualitative observations 

(1) As the oscillatory region moves to larger k, it "eats" Fermi surfaces in the (k, q) 
plane. These Fermi surfaces branches move to higher q for larger \g m \ (see Figure 




Fig. 2. Fermi Surfaces for Increased \g m \ with m = 0.4 



(2) The dipole coupling has the most effect at low q, where it flattens and curves 
Fermi surface branches close to the oscillatory region. Far from the oscillatory 
regions, the branches asymptote to straight lines. 

(3) This effect is most pronounced for m < (alternative quantization). For m 
negative enough, local and global maxima and minima can develop in Fermi 
surface branches near the oscillatory region. See Figure [3] for an example of how 
such a minimum develops as m is lowered. Also, in Figure [4] we plot results for 
g m fixed and m = —0.4, 0, 0.4. 



In 14 , a "phase diagram" was constructed in the (m, q) plane which showed the 
attainable v^s for the primary Fermi surface (that with the largest kf for a given 
q). Here, we construct a similar phase diagram for g m = —2. Because of various 
ambiguities that arise when there are multiple Fermi surfaces, we focus on a single 
q > branch. Because the branch gets flattened near the oscillatory region, there 
are more attainable u's than with g m = for the same range of m, q. Also, because 
of point 3 above, when global or local minima occur we must pick what we mean 
as the primary Fermi surface within a branch (note that this differs from g m = 
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Fig. 3. Developing local and global minimum for m < with g m = —2 




Fig. 4. Fixed g m = -2 for m = -0.4, 0, 0.4 



where the choice is made between different branches) . We choose that Fermi surface 
with the largest v a . Note, in this case, such a Fermi surface actually has smaller kf. 




-OA -^.2 0.0 0.2 04 -0.4 -Q2 0.0 0.2 0.4 



Fig. 5. "Phase Diagram" for g m = —2 



6.2. g e j£ 0, g m = 

For (/ e ^ 0, there is mixing between the spin components. We change bases in the 
uj = Dirac equation so that ( 68 ) is the near horizon limit. We then use two different 
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infalling boundary conditions, each corresponding to a distinct AdS2 dimension. We 
integrate this out to the boundary, change basis back to the original spi n b asis and 




Fig. 6. Increasing \g e \ 



(1) Fermi surface branches continue to jut out of an oscillatory region. As \g e \ is 
increased, the oscillatory regions kiss and move to the right (see figure [6]). As 
in the magnetic case, the oscillatory region "eats" Fermi surface branches as it 
moves to the right. 

(2) Fermi surface branches are created to the left (smaller |fc|) of the oscillatory 
region, which we will call the "interesting region". As one increases \g e \ more 
Fermi surfaces are created in this region. There can also be local maxima or 
minima created near the oscillatory region as in the magnetic dipole case. 

(3) For g e ^ the Fermi surfaces are much more gently sloping in the interesting 
region than for Fermi surface branches with similar g m (and all other constants 
comparable). For m = 0, the Fermi surface branches are nearly flat. 

(4) There seem to be small gaps between the Fermi surface branches at k = 0. This 
indicates local maxima or minima at k = 0. As one lowers \m\, the gaps become 
larger and larger, although for m < there always exists a branch with a large 
gap. 

(5) Since g e mixes the spin components, it is difficult to represent the landscape of 
possible values of v in the same way as previously. 

7. Discussion 

We have found that the holographic framework for Fermi surfaces is robust under 
the change of the magnetic and electric dipole moments of the bulk spinor field, in 
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the sense that the low-frequency Green's function remains of the form determined 



m 



14 . Turning on these couplings moves the Fermi surfaces around, and changes 
the scaling dimension of the emergent IR AdS2 symmetry. 

What is the meaning of the bulk dipole couplings in terms of properties of the 
fermion operator in the boundary field theory? In vacuum, these couplings do not 
affect the fermion two point function. They do, however, change the structure of the 
current-fcrmion-fermion three point function. The full calculation is complicated, 
but one simple characterization is the following. With g m = g e = 0, there are no 
terms in (^> a ^i3 J* 1 ) proportional to second rank Clifford algebra elements turn- 
ing on the dipole couplings creates such terms. It would be interesting to understand 
better the physical significance of this. 

It was previously found that the existence of Fermi surfaces in ordinary quanti- 
zation was correlated with the existence of the oscillatory region in k space where 
the AdS2 scaling dimension becomes imaginary. The oscillatory region occurs for 
values of the momentum such that there is Schwinger pair production in the AdS2 
region. A heuristic interpretation is that this creates a bulk Fermi surface, leading 
to a boundary one (the duality implies an equality of Hilbert spaces, meaning that 
the bulk spectral density is related to the boundary spectral density; sed 15 l 24 l for 
the precise relation between these two quantities). It appears that this mechanism 
involving pair production is not necessary to have a Fermi surface; boundary Fermi 
surfaces can appear - in alternative quantization - without oscillatory regions. How- 
ever, the alternative quantization is unstable in the RG sense; any small addition 
of the double trace operator O'O flows the CFT to that of ordinary quantization. 
Thus, the known Fermi surfaces without oscillatory region are also unstable; they 



flow away, as can be seen in Figure 5 of 14 



The existence of an oscillatory region implies that the bulk Fermi sea has sup- 
port at the black hole horizon. Naively, its gravitational backreaction should be 
suppressed by a factor of the Newton's constant, which is proportional to 1/N 2 . 
However, when one integrates the near-horizon charge density to some radial posi- 
tion r, there is a logarithmic divergence in r which can offset this suppression beyond 
some critical r c ~ e~ N ; for r < r c backreaction cannot be ignored.^ This backre- 
action was argued to change the geometry to a Lifshitz geometry^ with dynamical 
exponent z ~ ./V 2 . This modifies the fermion response for frequencies and tempera- 
tures below some new low-energy scale E c ~ fie~ N ~~- Below these frequencies and 
temperatures, the behavior is that of a Fermi liquid (the self-energy is analytic in 
frequency at the Fermi surface), and therefore the system does not in fact describe a 
non-Fermi liquid groundstate. We note that it is likely that other instabilities, such 
as the holographic superconductor instability^^ se t in at much higher tempera- 
tures. The effects of such superconducting order on holographic Fermi surfaces has 



been studied first for s-wave order 32 ■ 34 and more recently for order parameters 
with nodes 16|35 ■ 38 . 

For the reason described in the previous paragraph, it would be interesting to 
find a deformation of the action which allows Fermi surfaces that are RG stable 
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but do not come with an oscillatory region. Unfortunately, our results indicate that 
the dipole couplings are not such a deformation. Although the couplings change the 
shape and location of the oscillatory region in the (fc, q) plane, we again find that 
Fermi surfaces in ordinary quantization only occur with an oscillatory region. 

A subject of current interest is the form of the gravitational backreaction of 
the density of spinor particles in the bulk HUMES Given that the dipole couplings 
studied here change the shape of the Fermi surface and of the oscillatory region in 
momentum space, they will also have an interesting effect on the character of the 
back-reacted solution^] 

Finally, we have found that the dipole operators curve Fermi surface branches 
in the (fc, q) plane close to the oscillatory region. For certain values of g m and g e we 
can create local maxima and minima of these branches. It would be interesting if we 
could embed this system into one where q is a tunable parameter. In this context, 
a local maximum, for example, would represent two Fermi surfaces that merge and 
annihilate as q is continuously increased. We leave such an embedding to future 
work. 

Note Added 3: The papers |44|45| attempt to interpret the motion (upon in- 
creasing g m ) of the Fermi surface pole into the oscillatory region and its subsequent 
evolution as the formation of a garj^j At first glance, this interpretation is prob- 
lematic, because a characteristic feature of the oscillatory region is nonzero spectral 
weight at zero frequency (see figure 4 of [12] and section IV. A of [14]). At values of 
the dipole coupling when the Fermi surface first disappears, it enters the oscillatory 
region and the IR CFT exponent becomes i times a small number; the low-frequency 
spectral weight, while incoherent, is by no means small there. The interpretation 
of 44|45 of the formation of a gap relies on the numerical smallness of this weight 



at still-larger values of dipole coupling. 

A partial explanation of this effect is the following. As the dipole coupling is 
increased further, for fixed k, one exits the oscillatory region again, and the IR 
CFT exponent (62 1 v = y/—q 2 /2 + m 2 L 2 + (k ± Cdg m ) 2 becomes real and positive 



and, eventually, large. Such a large, real IR CFT scaling dimension v suppresses 



the incoherent spectral weight away from poles of the Green's function 33 . This is 



because the low frequency spectral density satisfies 14 



2v 

oj ' 



ImG fl oc(-) <1 (71) 



g We thank David Vegh for discussions of this issue. 

h JM would like to thank the authors of 44|45 for correspondence which led to the following more 
careful consideration of these points. 
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for uj < fx and v 3> 1. [][j][jThis explanation is only a partial one because for given 
g m , there is still an oscillatory region (at k oc g m ), where there could in principle 



be gapless excitations. 45 observes that even in that regime, where the IR CFT 



dimensions are i times an order-one number, the spectral weight is suppressed. An 

,(o) ,(0) 



understanding of this effect must involve the behavior of the UV coefficients a± , 



with g„ 



Despite this analytic partial understanding of the phenomenon observed in 44 



45 , we retain some reservations about the interesting proposed connection to Mott 
physics. In particular, the oscillatory region is playing a crucial role in destroying the 
Fermi surface in this discussion. Much about its interpretation remains mysterious, 
particularly in light of its implications^! for gravitational back-reaction. 

Finally, we note that the same situation of no Fermi surfaces and large IR CFT 
dimension can also be reached by increasing mL, the mass of the spinor field in units 
of the AdS radius, at fixed charge and no dipole coupling; in this regime (m > qea) 
there is no oscillatory region. 
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